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Abstract 

During the process of thermal leptogenesis temperature decreases by about one order 
of magnitude while the baryon asymmetry is generated. We present an analytical 
description of this process so that the dependence on the neutrino mass parameters 
becomes transparent. In the case of maximal CP asymmetry all decay and scattering 
rates in the plasma are determined by the mass Mi of the decaying heavy Major ana 
neutrino, the effective light neutrino mass fhi and the absolute mass scale fn of the 
light neutrinos. In the mass range suggested by neutrino oscillations, mgoi — 8 x 
10~^ eV < rhi < rriatm — 5 x 10~^ eV, leptogenesis is dominated just by decays and 
inverse decays. The effect of all other scattering processes lies within the theoretical 
uncertainty of present calculations. The final baryon asymmetry is dominantly 
produced at a temperature Tg which can be about one order of magnitude below 
the heavy neutrino mass Mi. We also derive an analytical expression for the upper 
bound on the light neutrino masses implied by successful leptogenesis. 



1 Introduction and summary 



Leptogcncsis [1] provides a simple and elegant explanation of the cosmological matter- 
antimatter asymmetry. A beautiful aspect of this mechanism is the connection between 
the baryon asymmetry and neutrino properties. In its simplest version leptogenesis is dom- 
inated by the CP violating interactions of the lightest of the heavy Majorana neutrinos, 
the seesaw partners of the ordinary neutrinos. The requirement of successful baryogenesis 
yields stringent constraints on the masses of light and heavy neutrinos. In particular, all 
light neutrino masses have to be smaller than 0.1 eV [2]. 

Leptogenesis is closely related to classical GUT baryogenesis [3], where the deviation 
of the distribution function of some heavy particle from its equilibrium distribution pro- 
vides the necessary departure from thermal equilibrium. The non-equilibrium process 
of baryogenesis is usually studied by means of Boltzmann equations [4, 5]. In the same 
way, leptogenesis has been studied during the past years, with increasing sophistication 
[6]- [12]. The goal of the present paper is to provide an analytical description of the lep- 
togenesis process such that the dependence on the neutrino mass parameters becomes 
transparent. As we shall see this is important to understand the size of corrections to the 
simplest Boltzmann equations, which have to be taken into account to arrive eventually 
at a 'theory of leptogenesis'. 

We will first consider the simplest case where the initial temperature Ti is larger than 
Ml, the mass of the lightest heavy neutrino A^i. We will also neglect decays of the two 
heavier neutrinos N2 and A^3, assuming that a generation oi B — L asymmetry from their 
decays either does not occur at all or that it does not influence the final value B — L. 
Further, we restrict ourselves to the non supersymmetric case, and we assume that the 
lightest heavy neutrino A^i is the only relevant degree of freedom beyond the standard 
model particle species. 

Within this minimal framework the Boltzmann equations can be written in the fol- 
lowing form^. 



dz 

dNB-L 



= -(D + 5)(Ar^,-Ar^^J, (1) 
= -eiD{Nj,,-N'j^^)-WNB-L, (2) 

where z = Mi/T. The number density Nj^^ and the amount oi B — L asymmetry, Nb-l, 
are calculated in a portion of comoving volume that contains one photon at temperatures 
T ^ Ml, so that the relativistic equilibrium Ni number density is given by N^^{z <C 
^We use the conventions of Ref. [10]. 
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1) = 3/4. There are four classes of processes which contribute to the different terms 
in the equations: decays, inverse decays, AL = 1 scatterings and AL = 2 processes 
mediated by heavy neutrinos. The first three all modify the A^i abundance and try to 
push it towards its equilibrium value iV^. Denoting by H the Hubble expansion rate, 
the term D = T£,/{H z) accounts for decays and inverse decays, while the scattering term 
S — Vs/{H z) represents the AL = 1 scatterings. Decays also yield the source term for 
the generation of the B — L asymmetry, the first term in Eq. (2), while all other processes 
contribute to the total washout term W = Tw/{H z) which competes with the decay 
source term. The expansion rate is given by 



where — qsm — 106.75 is the total number of degrees of freedom, and Mpi = 1.22 x 
10^^ GeV is the Planck mass. Note that we have not included the A^i degrees of freedom 
since, as we will see, in the preferred strong washout regime, the heavy neutrinos are 
non-relativistic when the baryon asymmetry is produced. 

The two terms D and S depend on the effective neutrino mass [7], defined as 

(4) 

Ml ^ ^ 

which has to be compared with the equilibrium neutrino mass 

= —— ~ 1.08 X 10"^ eV . 5 

3V5 Mpi ^ ^ 

The decay parameter 

_ Fpiz ^ oo) _ mi 
H{z = 1) m, ' 

introduced in the context of ordinary GUT baryogenesis [3] , controls whether or not A^i 
decays are in equilibrium. Here Tuiz = oo) = Fd is the A^i decay width. The washout 
term W has two contributions, W = Wq + AW; the first term only depends on rhi, while 
the second one depends on the product Mi m^, where = mf + m| + m| is the sum of 
the light neutrino masses squared [10]. 

The solution for Nb-l is the sum of two terms [3] , 

Nb-l{z) = N'b_^ e-^^ ''^'^(^'^ - ^ £i k{z; mi. Mi m') , (7) 

where the second term describes B — L production from A^i decays. It is expressed in 
terms of the efficiency factor k [9] which does not depend on the CP asymmetry ei. In 



the following sections we shall use two integral expressions for the efficiency factor, 

^{z) = ^ r dz'D{N^, - N^^J e-r, <^^"^^^") (8) 

~ 3 D + S dz' ■ ^ ' 

Here N^^ and dN^^/dz' are the solution of the first kinetic equation (1) and its derivative, 
respectively. The efficiency factor k{z) is normalized in such a way that its final value 
Ki = k{oo) approaches one in the limit of thermal initial abundance of the heavy neutrinos 
A^i and no washout (W = 0). In general, for N}^^ < N'^_^ =3/4, one has < 1. The first 
term in Eq. (7) accounts for the possible generation oi a. B — L asymmetry before Ni 
decays, e.g. from decays of the two heavier neutrinos N2 and A^3, or from a completely 
independent mechanism. In the following we shall neglect such an initial asymmetry 
N^B-L- [2] it was shown that for values rhi > m* even large initial asymmetries are 
washed out for initial temperatures Ti > Mi. 

The predicted baryon to photon number ratio has to be compared with the value 
TjB measured at recombination. It is related to iV^_^ = Nb-l{z = 00) by riB = 
(o-sph//) ^B-L- Here flgph = 28/79 [13] is the fraction oi B — L asymmetry converted 
into a baryon asymmetry by sphalcron processes, and / = N^^'^/N* = 2387/86 is the dilu- 
tion factor calculated assuming standard photon production from the onset of leptogenesis 
till recombination. Using Eq. (8), one then obtains 

rjs^ ^£^Ki = deiKtc^ 0.96 x 10~^ SiKf . (10) 
4 / 

In the following sections we will study analytically the solutions of the kinetic equa- 
tions, focusing in particular on the final value of the efficiency factor. We start in sect. 2 
with the basic framework of decays and inverse decays. In the two regimes of weak 
(mi < m*) and strong (mi > m*) washout the efficiency factor is obtained analytically, 
which then leads to a simple interpolation valid for all values of mi. AL — 1 scatterings 
are added in sect. 3, and the resulting lower bounds on the heavy neutrino mass Mi 
and on the initial temperature Ti are discussed. In sect. 4 an analytic derivation of the 
upper bound on the light neutrino masses is given, and in sect. 5 various corrections are 
described which have to be taken into account in a theory of leptogenesis. In appendix 
A a detailed discussion of the AL = 2 processes in the resonance region is presented. 
In the case of maximal CP violation the entire AL = 2 scattering cross section can be 
expressed in terms of Mi, mi and m. The resulting Boltzmann equations are compared 
with previously obtained results based on exact Kadanoff-Baym equations. In appendix 
B various useful formulae are collected. 
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Recently, two potentially important, and usually neglected, effects on leptogenesis 
have been discussed: the AL — 1 processes involving gauge bosons [11, 12] and thermal 
corrections at high temperature [12]. Further, the strength of the AL = 2 washout term 
has been corrected [12] compared to previous analysis. However, the reaction densities 
for the gauge boson processes are presently controversial [11, 12]. Also the suggestion 
made in [12] to include thermal masses as kinematical masses in decay and scattering 
processes leads to an unconventional picture at temperatures T > Mi, which differs 
qualitatively from the situation at temperatures T < Mi. If thermal corrections are 
only included as propagator effects [14] their influence is small. This issue remains to be 
clarified. Fortunately, both effects are only important in the case of weak washout, i.e. for 
mi < m^, where the final baryon asymmetry is strongly dependent on initial conditions 
in any case. In the strong washout regime, mi > m*, which appears to be favored by 
the present evidence for neutrino masses, they do not affect the final baryon asymmetry 
significantly. In the following we will therefore ignore gauge boson processes and thermal 
corrections. These questions will be addressed elsewhere. 

The main results of this paper are summarized in the figures 6, 9 and 10. Fig. 6 
illustrates that for the basic processes of decays and inverse decays the analytical approx- 
imation for the efficiency factor agrees well with the numerical result. The figure also 
demonstrates that scatterings lead to a departure from this basic picture only for val- 
ues K — mi/m^ < 1, where the final baryon asymmetry depends strongly on the initial 
conditions. Fig. 9 shows the dependence of the efficiency factor on initial conditions and 
on AL = 1 scatterings for different values of the effective Higgs mass M^. Again, for 
mi > m* this dependence is small and, within the theoretical uncertainties, the efficiency 
factor is given by the simple power law 



Knowing the efficiency factor, one obtains from Eqs. (10) and (118) the maximal baryon 
asymmetry. Fig. 10 shows the lower bound on the initial temperature T; as function 
of mi. In the most interesting mass range favored by neutrino oscillations it is about 
one order of magnitude smaller than the lower bound on Mi. The smallest temperature 
jjmm ~ 3 X 10^ GeV is reached at mi ~ 2 x 10~^ eV. In sect. 4 an analytic expression for 
the light neutrino mass bound is derived, which explicitly shows the dependence on the 
involved parameters. 




(11) 
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2 Decays and inverse decays 



It is very instructive to consider first a simplified picture in which decays and inverse de- 
cays are the only processes. For consistency, also the real intermediate state contribution 
to the 2 — > 2 processes has to be included. The kinetic equations (1) and (2) then reduce 
to 



dz 
dNs-L 
dz 



= -^(TV^.-TV^^), (12) 
= -SiD {Nm,-N''^J-WjdNb-l, (13) 



where Wm is the contribution to the washout term due to inverse decays. Prom Eqs. (8) 
and (12) one obtains for the efficiency factor, 

^(^) ^ _1 r dz' ^ e-r.' dz"w,u{z") (14) 
3 dz' 

As we shall see, decays and inverse decays are sufficient to describe qualitatively many 
properties of the full problem. 

After a discussion of several useful analytic approximations we will study in detail 
the two regimes of weak and strong washout. The insight into the dynamics of the non- 
equilibrium process gained from the investigation of these limiting cases will then allow 
us to obtain analytic interpolation formulae which describe rather accurately the entire 
parameter range. All results will be compared with numerical solutions of the kinetic 
equations. 

2.1 Analytic approximations 

Let us first recall some basic definitions and formulae. The decay rate takes the form [4], 

Toiz) = , (15) 

where the thermally averaged dilation factor is given by the ratio of the modified Bessel 
functions Ki and K21 

1\ K.i^) (16) 



.7/ M^) 

and To is the decay width. 



- "8^^ ' ^^^^ 
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Figure 1: The dilation factor. The dashed line is the analytical expression Eq. (26) to be 
compared with the numerical result (solid line). 



with the Higgs vacuum expectation value v = 174 GeV. The decay term D is conveniently 
written in the form [3] 



D(z) =Kz 
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The inverse decay rate is related to the decay rate by 

Tio{z) = Td{z) ' 



(18) 



(19) 



where is the equilibrium density of lepton doublets. Since the number of degrees of 
freedom for heavy Majorana neutrinos and lepton doublets is the same, g^-^ = gi = 2, one 
has 

3 2 



The contribution of inverse decays to the washout term W is therefore 

1 Tjd(z] 



Winiz) 



(20) 



(21) 



which, together with Eqs. (16), (18) and (20), implies 



Wjniz) = -D{z) 



(22) 
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All relevant quantities are given in terms of the Bessel functions Ki and K2, whose 
asymptotic limits are well known. At high temperatures one has, 

K2(z) ~ - K,(z) ~ 4 ' < 1 ' (23) 
z z 

whereas at low temperatures, 

Accurate interpolating functions for Ki{z) and K2{z) for all values of z are 



Note, that for 2; <S 1 this approximation gives K2{z) ~ 15/(82;^) rather than the exact 
asymptotic form 2/z'^. However, the high temperature domain is not so important for 
baryogenesis and the approximation (25) is rather precise in the more relevant regime 
around 2; ~ 1. 

Eq. (25) yields very simple expressions for the dilation factor and the decay term, 

DU)^K^^. (26) 

As Fig. 1 shows these analytical approximations are rather precise. The relative error is 
always less than 7%. The washout term (21) becomes in the approximation (25), 



Wjn{z)^lKz'^l + ^ze-^ . (27) 

It is useful to define a value z^, corresponding to a 'decay temperature' below which 
decays are in equilibrium, by rD{zd)/H{zd) — ZdD{zd) — 2. The value of z^ is determined 
by K, and from Eq. (26) one obtains 



For X < 1, this yields za ^ •\/2/X, whereas z^ ~ (15/4ir)^/^ for X > 1. At X ~ 1 one 
has Zd ~ 2. 

Inverse decays are in equilibrium if Wid{z) > 1. From Eq. (27) one easily finds that 
Wid{z) reaches its maximal value Wioizraax) — 0-SK at Zmax — 2.4. Hence, for K > 3, 
there exists an interval < z^^x < Zout, where inverse decays are in equilibrium. For 
K <3 no such interval exists and inverse decays are always out of equilibrium. 




Figure 2: Out of equilibrium decays. Ni number density, efficiency factor and decay 
temperature = Mi/z^ for K = 10~^ and K = 10~^. 

2.2 Out-of-equilibrium decays 

In the regime far out of equilibrium, K <^ 1, decays occur at very small temperatures, 
^ 1, and the produced B ~ L asymmetry is not reduced by washout effects. In this 
case the integral for the efficiency factor (14) becomes simply, 

• (29) 

For 2 < no asymmetry is generated because the heavy neutrinos do not decay. 
They also cannot be produced since inverse decays are switched off as well. Hence, in 
this regime the dynamics is completely frozen. For z > za the equilibrium abundance is 
negligible, and from Eq. (12) one finds, 

^ iVi,^e-^(--^^(¥)^^"(^-^*^)). (30) 

Note that we have neglected the small neutrino abundance which for N\j^ <^ is 
produced before the neutrinos decay. Fig. 2 shows the evolution of N^^i^z) and Nb-l{.z) 
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for K — 10~^ and K — 10~^, with N]^_^ — =3/4, comparing the numerical solutions 
with the analytical expressions. 

The final value of the efficiency factor Kf — k{oo) is proportional to the initial A^i 
abundance. If = A^^'^ = 3/4, then Kf = 1. But if the initial abundance is zero, then 
= as well. Therefore in this region there is the well known problem that one has to 
invoke some external mechanism to produce the initial abundance of neutrinos. Moreover 
the assumption that the initial asymmetry is washed out does not hold. Thus in the 
regime K <^ 1 the results strongly depend on the initial conditions and the picture is not 
self-contained. 



2.3 Dynamical initial abundance 

In order to obtain the efficiency factor in the case of vanishing initial A^i-abundance, 
A^jVi (zi) = -/Vjvi — 0> calculate how heavy neutrinos are dynamically produced 

by inverse decays. This requires solving the kinetic equation (12) with the initial condition 
Nh, = 0. 

Let us define a value z^q by the condition 

NN^^e^) = K\{Z^) . (31) 

Eq. (12) implies that the number density reaches its maximum aX z — Zeq- An approximate 
solution can be found by noting that for z < Zeq inverse decays dominate and thus 



dz ^^^^^1 



~ L> A^^^ > . (32) 



A straightforward integration yields for z < Zeq (cf. (16), (18), (20)), 

N^,{z) ~ £dz'z"K,{z') 

= \£ dz'Wioiz'). (33) 
In the case zi <^ z < 1, this implies for the number density, 

Nr,,{z)^jz', (34) 

where the small dependence on has been neglected. 

We can now calculate the corresponding approximate solution for the efficiency factor 
k{z). For z < Zeq the efficiency factor «; = ac~ is negative since, N^-^ < . From 
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Eqs. (14) and (32) one obtains 

K- {z) ~ ~ £ dz'D{z')N^^^ {z') e- <^-"^'o{z") 

= -2 r dz'Win{z')e- ^^'"'""^^-^'"^ 

~ -2(1- e-i , ^ < ^eq ■ (35) 

As expected, for A^jvi <^ 1 the efficiency factor is proportional to A^at^, up to corrections 
whicfi correspond to washout effects. For z > Zeq, k'~{z) is reduced by washout effects. 
For z > Zeq, there is an additional positive contribution to the efficiency factor, 

K+{z) ~ - r dz'D{z') {NnAz') - N'lliz')) e-^'"^''"^'^^'"^ . (36) 

The total efficiency factor is the sum of both contributions, 

Kl{z) = {z) ^ K- {z) . (37) 

For z > Zeq we now have to distinguish two different situations, the weak and strong 
washout regimes, respectively. 

2.3.1 Weak washout regime 

Consider first the case of weak washout, -C 1, which implies z^q I- From Eq. (33) 
one then finds, 

Qtt 

N^.izeq)- — K = N{K) . (38) 
16 

A solution for Nisf^{z), valid for any z, is obtained by using in Eq. (33) the useful approx- 
imation 

37r 7^ 

dz' z'^ KAz')^ , , (39) 

^ ' [(97r)'=+ (2^3)c]i/c ' ^ > 

with c = 0.7. This yields an interpolation of the two asymptotic regimes (cf. Eqs. (34) 
and (38)), which is in excellent agreement with the numerical result, as shown in Fig. 3a 
for K = 10-2. 

For 2; > ^eq ^ 1 decays dominate over inverse decays, such that 



dATjVi 
dz 



-D Nn,<0. (40) 
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In this way one easily obtains 



A^iv.(^)=A^^,We-^i.'^^'^(^r (41) 

Moreover, for z > z^^, Wid{z) is exponentially suppressed and washout effects can be 
neglected in first approximation. For the negative part of the efficiency factor one then 
has (cf. Eq. (35)) 

K'iz) = -2 (l - e-i^(^)) . (42) 
Prom Eq. (29) one obtains for the positive contribution, 

k+{z)^^{N{K)-Nn,{z)) . (43) 

The final efficiency factor is then given by 

Kf{K) ~ ^ N{K) -2(1- e-i^(^)) . (44) 

To first order in N{K) oc K the final efficiency factor vanishes. This corresponds to 
the approximation where washout effects arc completely neglected. As discussed above, 
Kf is then proportional to iVj^^ and therefore zero. To obtain a non-zero asymmetry the 
washout in the period z < Zeq is crucial. It reduces the absolute value of the negative 
contribution k,~{z), yielding a positive efficiency factor of order 0{K'^), 



\n{k) 



^^X^ (45) 



Such a reduction of the generated asymmetry has previously been observed in the context 
of GUT baryogenesis [15]. Note, that for X > 1 Eq. (45) does not hold, since in this case 
2;eq becomes small and washout effects for 2; > are also important. 

In Fig. 3a the analytical solutions for Nn^{z) and \Nb-l{z) \ — (3/4)|£i«;(2;)| are com- 
pared with the numerical results for K ~ 10~^. A residual asymmetry survives after 
z^ ^ 2;eq as remnant of the cancellation between the negative and the positive contribu- 
tions to the efficiency factor. The second one is prevalent because washout suppresses 
K~ more efficiently. As one can see in Fig. 3a, the analytical solution for the asymme- 
try slightly underestimates the final numerical value. This is because for K > 10~^ the 
approximation of neglecting washout for z > Zeq becomes inaccurate. 

2.3.2 Strong washout regime 

As K increases, Zeq decreases, and at X ~ 3 the maximal number density N{K) reaches 
the equilibrium density Neq at Zeq ^ 1. For K ^ 1, one obtains from Eq. (34), Zeq ~ 
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Figure 3: Comparison analytical (dashed lines) and numerical (solid lines) results for 
heavy neutrino production and B — L asymmetry in the case of zero initial abundance, 
= 0, for weak washout (top) and strong washout (bottom); = 10~^. 
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(e/X)^/^ ~ <^ 1- A more accurate description for z < Zeq has to take into account 
decays in addition to inverse decays, i.e. one has to solve Eq. (12). Since -Zeq ^ 1, one 
can use — 3/4, and one then easily finds, 

N^,{z)^\(\-e-\^^'), (46) 

which correctly reproduces Eq. (34) for z <^ 2;eq. An example, with K — 100, is shown in 
Fig. 3b which illustrates how well the analytical expression for neutrino production agrees 
with the numerical result. 

Consider now the efficiency factor. For X ^ 1 we can neglect the negative contribution 
K", assuming that the asymmetry generated at high temperatures is efficiently washed 
out. This is practically equivalent to assuming thermal initial abundance. We will see 
in the next section how to describe the transition from the weak to the strong washout 
regime. 

For K > 3, inverse decays are in equilibrium in the range z\a < z < z^^x^ with 
2;in 2/\fK. In the strong washout regime, K ^ 1, the efficiency factor can again be 
calculated analytically. 

For z < za — [15/(4 X)]^/^ decays are not effective in tracking the equilibrium distri- 
bution. For the difference 

A = 7V^,(^)-7V^^(^), (47) 
with A^jVi(O) = A^^j(O) = 3/4 = A^eq, one has, 

/^{z)^N'^^-N^^^{z)^^^z\ z<z,. (48) 

The corresponding efficiency factor is given by (cf. (8)) 

4 2K 
K{z)^-j dz'D{z')A{z')^—z\ z<za. (49) 

On the other hand, for z > z^ the neutrino abundance tracks closely the equilibrium 
behavior. Since D oc K, one can solve Eq. (12) systematically in powers of l/K, which 
yields 

Using the properties of Bessel functions, Eq. (20) yields for the derivative of the equihb- 
rium density. 
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We can now calculate the efficiency factor. Prom Eqs. (14) and (51) one obtains 

K J^. z' 

dz' e-"^^''''^ . (52) 



Zi 



The integral is dominated by the contribution from a region around the value z where 
ip{z', z) has a minimum. The condition for a local minimum zb, the vanishing of the first 
derivative, yields 

Wjuiz^) = (-) {zb) - - . (53) 

\7/ 

Since the second derivative at zb is positive one has z = minj^;, zb}- 

The integral (52) can be evaluated systematically by the steepest descent method 
(cf. [3]). Alternatively, a simple and very useful approximate analytical solution can be 
obtained by replacing in the exponent of the integrand Wid{z) by 

Wid{z) = ^ Wjd{z) = ^(^'^2(;^)) . (54) 
The efficiency factor then becomes 

«;(;,) ~ ^ r dz'WiD{z')e-^^^"''"^'^^'"^ 

A fi_e-/4'^^'^-(^'))) . (55) 
Kz V / 



It is now easy to understand the behavior of k{z). For z^ < z < zb, one has k oc l/z, 
while for z > zb the efficiency factor gets frozen at a final value Ki ~ 2/{K zb), up to a 
small correction 0(exp [—K)). 

One can also easily find global solutions for all values of z by interpolating the asymp- 
totic solutions for z < z^^^ and z > z^, respectively. From Eqs. (18), (48), (50) and (51) 
one obtains for the difference between A^i-abundance and equilibrium abundance, 

A(^)-(l + y:^j j^z'Mz). (56) 

Similarly, an interpolation between the expressions (49) and (55) for the efficiency factor 
is given by 

, , / K''zz^Y^2K 5 
-(^)-(l + ^j 7^^- (^^) 

A typical example of strong washout is shown in Fig. 4 for the value K — 100, as in 
Fig. 3b, but now for thermal initial abundance. In this figure we also show the decay, 
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Figure 4: Strong washout: comparison between analytical and numerical (full lines) results. 
Inverse decays are in equilibrium in the temperature range ^in < < -^out ~ = 10~^. 

inverse decay and washout terms. Instead of the neutrino abundance the deviation A (2;) 
is shown. The dotted, short dashed, dot-dashed and dashed lines are the approximations 
Eq. (26) for the D term, Eq. (27) for the Wiy) term, Eq. (56) for A(z) and Eq. (57) for 
fi:(z), respectively. The thin solid lines are the numerical results which agree well with 
the analytical approximations. The behavior k{z) (xl/ z for z > and the freeze-out of 
Nb-l at Zo^t are clearly visible. 

Let us now focus on the final value of the efficiency factor Kf = n{oo). Note, that for 
K ^ 1 also z-Q ^ 1, and the condition (53) becomes approximately W/d(-2b) — 1- This 
means zb — Zont- Hence, the asymmetry produced for z < Zout is essentially washed out, 
while for for z > Zou_t washout is negligible {Wid < 1). This simple picture will have some 
interesting consequences and applications. 

The integral in Eq. (55) is easily evaluated, 

dz'WiDiz) = ^ZBiK)K . (58) 
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Figure 5: zb as function of the decay parameter [3K. The case studied in this section 
corresponds to /5 = 1. The sohd red hne is the numerical solution of eq. (59), the dotted 
black line shows the asymptotic solution (60), and the blue dashed line is the interpolation 
(62). 

Using the approximations (18) and (27), the condition (53) for z-q{K) becomes explicitly 
^ZB(i^)'e-^B(^) ^1 + '^z^{K) ^ z^{K) - 1 . (59) 
Zb{K) approaches one as K goes to zero^. For K ^ 1 the solution of Eq. (59) is given by 



57rV5 



(60) 



where W-i is one of the real branches of the Lambert W function [16]. This result can 
be approximated by using the asymptotic expansion of W^i [16, 17], 



2 r 



ttK 

Zb(K) ^ - In I 

^ ^ 2 I 1024 



In 



/31257rJs:2 
V 1024 



+ 



\n{K) 
ln(ln(ir)) 



(61) 



A rather accurate expression for zb{K) for all values of K is given by the interpolation 
(cf. Fig. 5), 



ZBiK)c^l + -ln I 1 



1 . I . nK^ r. /31257rfs:2x n5 



1024 



In 



V 1024 



(62) 



^Note that the solution zb(K) of Eq. (53) approaches asymptotically 1.33 for K 0. However, in 
the strong washout regime and also for the extrapolation K this difference is irrelevant for Kf. 
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Figure 6: The final efficiency factor Kf as function of tlie decay parameter K for thermal 
and dynamical initial A^i abundance, respectively. The thick solid lines are the numerical 
solutions. The thin lines show, for comparison, the numerical solutions of the complete 
Boltzmann equations including A'^i-top scatterings with an effective Higgs mass M^/Mi = 
0.1. The red circles represent Eq. (63) for the efficiency factor, evaluated using Eq. (62) 
for zsiK). 

Note the rapid transition from strong to weak washout at ^ 3. 
The final value of the efficiency factor takes the simple form 



This analytical expression for the final efficiency factor, combined with Eq. (62) for 
zb{K), provides an accurate description of K,f{K), as shown in Fig. 6. Eq. (63) can also 
be extrapolated into the regime of weak washout, K <^ 1, where one obtains Hf = 1 
corresponding to thermal initial abundance, N}^_^ = = 3/4. It turns out that also in 
the transition region Eqs. (63) and (62) provide a rather accurate description, as is evident 
from Fig. 6. The largest discrepancy between analytical and numerical results is about 
30% around if ~ 1. For comparison also the numerical result including scatterings is 
shown. The difference with respect to the basic 'decay-plus-inverse-decay' picture becomes 
significant only for K < 1. 




(63) 
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The above analysis is easily extended to the case where the strength of the washout 
term Wm is modified to pWjD- For instance, in the model considered in [3], S — L number 
changes by two in heavy particle decays, corresponding to /? = 2. On the other hand, 
the heavy particle abundance is not afi^ected by this change. The final efficiency factor is 
therefore given by 

- --™) • 

where Zb{K) is again given by Eq. (59). In the regime ^ 1 our expression for the final 
efficiency factor can be approximated by 

2 1 

'^'^-^^ - zb{K)(5K - 1.2 (5K {hi (5Kf-^ ' ^^^^ 

which is very similar to the result^ obtained by Kolb and Turner [3] . 

Comparing the efficiency factor (63) with the solution of the Boltzmann equations 
including scatterings, as shown in Fig. 6, one arrives at the conclusion that the simple 
decay-plus-inverse-decay picture represents a very good approximation for leptogenesis in 
the strong washout regime. As we will see, the difference is essentially negligible within 
the current theoretical uncertainties. 

2.3.3 Global parametrization 

Given the results of the previous sections it is straightforward to obtain an expression for 
the efficiency factor for all values of K also in the case of dynamical initial abundance. 
We first introduce a number density N{K) which interpolates between the maximal num- 
ber densities A^eq = 3/4 and N{K) — dnK/lS (cf. (38)) for strong and weak washout, 
respectively, 

= ^ . (66) 

The efficiency factor is in general the sum of a positive and a negative contribution. 

Here k^{K.) is given by (42) for -C 1. A generalization, accounting for washout also for 
z > Zeq, reads 

,lNiK)-J^' dz'Wiuiz') ^QY) 



K {Z} = K {Zcq ) e 



■^Quantitatively, a discrepancy by a factor ~ 7 was noted in [18], which is mostly related to the 
definition of the decay parameter K (cf. Eq. (6)). 
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This expression extends the vahdity of the analytical solution to values K > 10~^ in the 
case of a dynamical initial abundance. Kf is exponentially suppressed for X ^ 1. An 
interpolation, satisfying the asymptotic behaviors at small and large K, is given by 

KTiK) = -2 e-t^(^) (et^(^) - l) . (68) 

On the other hand, the expression (63) for k,^{K), which is valid for K ^ 1, has to 
approach 4/3 N{K) for K <^1. These requirements are fulfilled by 

4iK) = (i _ e-§-(^)^^W) . (69) 

Equations (66), (68) and (69), together with the interpolation (62) for zb{K), yield an 
accurate description of the efficiency factor K,i{K) for all values of K, as demonstrated by 
Fig. 6. 

This result is a good starting point for obtaining an analytic description of the effi- 
ciency factor for the full problem. In the following sections we shall go beyond the simple 
decay-and-inverse-decay picture and include other processes step by step. 



3 The scattering term 
3.1 Analytic approximations 

The scattering term S and the related washout contribution Wal=i arise from two differ- 
ent classes of Higgs and lepton mediated inelastic scatterings involving the top quark (t) 
and gauge bosons (^4), 

S^S^ + S^ . (70) 

Their main effect is to enhance the neutrino production and thus the efficiency factor 
for mi < m*. Further, they also contribute to the washout term, which leads to a 
correction of the efficiency factor for rhi > m*, i.e. in the strong washout regime. Since 
the scattering processes are specific to leptogenesis we shall use in this section mostly 
the variable fhi > mi [19] rather than K. Top quark and gauge boson scattering terms 
are expected to be of similar size. However, the reaction densities for the gauge boson 
processes are presently controversial [11, 12]. We shall therefore discuss these processes in 
detail elsewhere. We shall also neglect the scale dependence of the top- Yukawa coupling, 
which reduces the size of S"*, since this decrease of S will be partially compensated by S^. 

The term is again the sum of two terms arising from the s-channel processes A^i I ^ 
t q and the t- channel processes Nit ^ Iq, Niq It, 

5* = 2V + 45<^,t. (71) 
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The scattering terms are defined as usual in terms of scattering rates and expansion rate, 



r(A^i) 
^ <t>At) 

Hz 



and introducing tlie functions /</,,s(t) (2;) (cf. appendix B) it is possible to write 



here we have introduced the ratio 



with 



mt ^ 



mi 
mt 



(72) 

(73) 
(74) 



9 mf 



10 m* . 



(75) 



At high temperatures, 2; <^ 1, the functions /0,s(t) have the following asymptotic form, 

'2 



l + iln( ^ 



-IE 



Ml 



In ( - ) - 7E 



(76) 
(77) 



In Fig. 7 wc have plotted the rates S + D and D as function of z for Ks = 1, i.e. 
"^1 = T^t- For values z < 2 the sum S* + -D is dominated by the scattering rate S, while 
for z > 2 the decay rate D c:^ K z dominates. A simple analytic approximation for the 
sum D + S is given by 



D + Sc^Ks 



l + ln(^UMn(l + - 



where 



K 



(78) 



(79) 



Ks HMjMh) 9 ln(Mi/M,) 
Here wc have introduced the Higgs mass Mh to cut off the infrared divergence of the 
t-channel process. As Fig. 7 illustrates, the approximation (78) agrees well with the 
numerical result for Mh/Mi = 10^^ as well as MfjMi = 10^^. Note that the latter value 
corresponds to the thermal Higgs mass, Mh — 0.4 T, at the baryogenesis temperature Tb 
in the strong washout regime. 

The washout term induced by the AL = 1 scattering processes is again the sum of s- 
and t-channel contributions. 



WAL=i^W^,, + 2W^,t 
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(80) 




z=M/T 



Figure 7: The rates S + D and D are shown as function of z for Ks = 1, i.e. mi = ml. The 
dash-dotted hue is the simple approximation S" + D ~ Ks + K z, while the two short dashed 
lines represent Eq. (78) for Mh/Mi = 10~^ (higher) and Mh/Mi = 10~^ (lower), to be compared 
with the numerical results (thick solid lines). The dotted line shows the expression (18) for D, 
to be compared with the numerical result for D (solid line). 

The washout rates are directly related to the scattering rates (72), 

^^,s - N'^^ H z ~ N'^\ ' ^ ' 

Using Eq. (22) for the equilibrium number densities one obtains 

Wal=i = 2WjD^ 5^,, + 2S^,t^ . (83) 

The washout rate including inverse decays is then given by 



Wo = Wid + Wal=i 
1 



W, 



ID 



1 + 



D 



M) 
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This is the total washout rate as long as AW, the off-shell contribution from heavy 
neutrinos, can be neglected. This is justified for sufficiently small values of Mi (cf. sect. 4). 



3.2 Dynamical initial abundance 

We can now calculate the production of heavy neutrinos and study how the efficiency 
factor gets enhanced by the presence of the scattering term. We again define a value Zeq 
by the condition (31), 

For z < Zeq the number density can be obtained by integrating the equation 

{D + S)N'^^ > . (85) 



The result is given by the expression {z < Zeq) 

- In ( 



N^,{z)^^Ks 



lJz)+ln(^]lB{z) 



(86) 



Here the first integral is given by 



37r 7^ 

lAz) = / dz' z- K,{z') ^ ^(g^^,^^^^3)c]Vc + ^^(-) ' (87) 







where we have used the approximation Eq. (39). The second integral can be expressed as 
Xb{z) = j'^ dz' z'^ In (l + ^) K2{z') 

~ 2 jT' dz' z'^ In (^1 + ^) + a dz' z"" K^iz') 

~ '^(^{l + a^)ln{l + a) -a^ In a- + ^a^ + a K3{1) -a z^ K3{z) . (88) 

The value Zeq can now be determined by setting A^jvi, as determined from Eqs. (86), (87) 
and (88), equal to N^^. Using an approximate form for , one obtains an equation 
similar to Eq. (59), as described in appendix B. This yields a good approximation for z^q 
in the case rhi < m*. 

3.2.1 Weak washout regime 

Consider now the case of weak washout, mi <C m* ~ 10~^ eV, which imphes 2;eq ^ 1. 
For z > Zeq, decays dominate over inverse decays. 



dNN, ^ 
dz 
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{D + S)Nn,<0. (89) 



Using D + 5" ~ Kz, valid for z ^ a (cf. (78) and fig. 7), this yields for the number 
density the simple expression 

Ar^^(^eq)e-^(^'-<). (90) 

In Fig. 8 the solution Nn-^(z) is shown for fhi — 10~^eV. The analytical solution 
agrees well with the numerical result. We also make a comparison with the result already 
displayed in Fig. 3, where the S term is neglected. As expected, the presence of the S 
term enhances the density A^^Vi at 2; = z^q- Moreover the comparison illustrates the strong 
sensitivity of the efficiency factor in the case mi <^ m*, not only to the initial conditions, 
but also to the theoretical description. A difference in A^i abundance by less than a factor 
of two at Zeq corresponds to final efficiency factors which differ by two orders of magnitude. 
This is due to delicate cancellations between the positive and the negative contribution 
to the efficiency factor and is a source of large theoretical uncertainties in the small rhi 
regime. 

We can now calculate the efficiency factor. In sect. 2.3.1 we have seen that in the 
absence of scatterings the inclusion of the small washout term was necessary to create 
an asymmetry between the negative contribution, k~, and the positive one, k'^, in order 
to have a non-zero final value Kf. Given the S term one can neglect washout to first 
approximation. From Eq. (8) one then obtains 

where (cf. (18), (78)) 



■( ^ ^ + ^ 



ilnfl + «U''^ 



a \ zJ Kz 



1 + ^1 . m 



Due to the S term we now have 7^ I^atj, although washout is neglected. The reason 
is quite clear: as in the case without scatterings, the asymmetry is changed only by 
decays and inverse decays; however, the number of decaying neutrinos at Zeq^ is now larger 
because of the additional production due to scatterings. To first approximation we can 
thus calculate the efficiency factor neglecting washout. 
For z < Zeq one obtains (cf. (39)), 

= -i £ i.'DN2 ^ -^^^^f^. . (93) 
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Figure 8: Effect of scatterings on neutrino production for fhi = 10^^ eV. The numerical results 
have been calculated for AW = and M^/Mi = 10~^ (thick solid lines). The short-dashed lines 
are the analytical solutions for A^at^ (Eqs. (86) and (90)) and for k{z) (Eqs. (93) and (94)), with 
Mh/Mi = 10"^ For comparison we also show the result where scatterings are neglected (thin 
solid lines). 



For z > Zeq one has j ~ const, (cf. Fig. 7); Eq. (91) the yields the simple result, 

k{z) = ^ (Ar^,(2:eq) J-^(Zeq) - N^A^) r\z)) + /^-(Zeq) , (94) 

which is shown in Fig. 8 for fhi = 10^^ and M^/Mi = 10~^ (short dashed line); it agrees 
reasonably well with the corresponding numerical solution (solid line). The analytical so- 
lution somewhat overestimates |; correspondingly, the final value k/ is underestimated. 

The analytical solution explains why the final value of the efficiency factor, Kf, is 
proportional to mi, 

4 



where 



iV^l(^eq) J-'(^eq) " N (K) OC . (95) 



2N{K)z, 



N{K) = , , (96) 



((97r)^ + (2^3^)^) 
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Figure 9: The final efficiency factor wiien scatterings are included. The numerical results 
are shown for M^/Mi = 10"^'', 10~^, 10~^, 1 (dot-dot-dash, solid, dashed and dot-dashed line 
respectively) . The thin dashed line is the simple result from the decay-plus-inverse-decay picture 
when scatterings are neglected. The short-dashed lines are the analytical results in the case of 
thermal initial abundance (Eq. (107) with j{zB) j'^{zB) in the exponential) and zero initial 
abundance (Eqs. (110) + (109), Mh/Mi = 10"^). The circled line is the power law fit (108). 
The dashed box indicates the range (msoi,w.atm)- 

with N{K) = dnK/W. Contrary to the case discussed in sect. 2.3.1, for which A^Ar^(zeq) — 
N{K) ~ N{K) and j = 1, NNi{zeq) and N{K) are now different. Hence, there is no 
cancellation of terms (9 (mi) between k,~^ and k,~. 

The expression (95) for the final efficiency factor fails for effective neutrino masses 
rhi > 10~^eV. Including washout mainly reduces the negative contribution k~ and 
thereby enhances the final value of the efficiency factor. Eq. (93) is then changed into 

^^i = -\ j dz' DN'^^^ e- . (97) 

For z < 2eq onc has NN-^ < A^^, and the washout rate becomes (cf. (84)) 

Wo{z)c:^Wjd{z)^^^^^ . (98) 
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Prom the expression (78) ior D + S one obtains 

..4.,,. ...(?.,n(^).M.(..3) 

~ '^K, + Kz; (99) 
here the last approximation requires a <S 1. Together with Eq. (18) this yields 

Wo{z)c^Wid{z)(i + ^) , (100) 
where the coefficient a is given by 

Since Wid/z is a total derivative, one obtains for the efficiency factor, 

POO 

POO 

^ -2 e-i^" / dz' Winiz') e"/^^ dz"w,oiz") ^ ^^32) 

where we have used z'^K2{z) ~ 2 for 2; < 1. Except for the exponential pre- factor, this 
yields the result obtained in sect. 2.3.1 for decays and inverse decays (cf. (42)) when the 
wash-out at 2; > z^ is neglected 

«;- = -2e-^^" (^l_e-§^(^)) . (103) 

The case without scatterings is recovered for a = 0. 

3.2.2 Strong washout regime 

In the case of strong washout, K ^ 1, the density of heavy neutrinos follows closely the 
equilibrium abundance, and one can obtain an analytical solution repeating the discussion 
in sect. 2.3.2. The efficiency factor is now given by 

4 r , , dN2 



K[Z) = — - 

^ ' 3 



r dz'3-' -^e-S^"'^"^°^^"K (104) 

J 2in dz 

Using iVjVi/iV^ ~ 1 in the washout rate one obtains (cf. (84)), 

W^^WiDj- (105) 
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In this way one finds for the efficiency factor 



(106) 



As in sect. 2.3.2, the dominant contribution to the integral arises from a region around 
a value zb ^ 1 where ips{z', z) has a minimum. Since -D + S* ~ Ks + Kz for large z, the 
value Zb is again given by Eq. (53) up to corrections 0{Ksl {K z^)). Replacing now Wid{z) 
by WiD{.z)zBi{zB) I {zj{z)) in the exponent of the integrand and by WiD{.z)j{z)/j{zB) in 
the pre-factor, respectively, one obtains for the final efficiency factor the approximate 
solution. 



This extends Eq. (64) to the case where scatterings are included. 

Note, that at small K the efficiency factor (107) does not approach one, the value corre- 
sponding to thermal initial abundance. However, any initial abundance can be reproduced 
by adjusting the exponent in Eq. (107). Replacing j{zB) by 3{zbY leaves Eq. (107) essen- 
tially unchanged at large K, whereas at small K one has ~ 1 corresponding to thermal 
initial abundance. The result is shown in Fig. 9 (short-dashed line) and compared with 
numerical results for different values of Mh/Mi. In the strong washout regime, mi ^ m*, 
and for Mh/Mi — 10~^, the analytical and numerical results agree within 10%. Since the 
strong washout regime is most interesting with respect to neutrino mass models, this is 
one of the most relevant results of this paper. 

For practical purposes it is interesting to note that, within the current theoretical 
uncertainties, the efficiency factor for rhi > m* is given by the simple power law. 



The quoted uncertainties represent, approximately, the range visible in Fig. 9 for large 
rhi. It is limited from above by the thin solid line corresponding to decays plus inverse 
decays and from below by the dot-dot-dashed line where scatterings are included with the 
extremely small ratio M^/Mi = 10~^°. Note that M^/Mi = 0.1 corresponds to a thermal 
Higgs mass, ~ 0.4 T, at the baryogenesis temperature Tb c± Mi/5. 

3.2.3 Global parametrization 

As in sect. 2.3.3 we can now obtain an expression for the efficiency factor for all values 
of K by interpolating between the two regimes of small K and large K. We shall use the 




(107) 




(108) 
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number density N{K) (cf. (96)) and the interpolation (62) for zb{K), which is related to 
the baryogenesis temperature by Tb — M-i/zb- 

The efficiency factor is the sum of a positive and a negative contribution, 

Ki{K) = Ki{K) + K^{K) . 

Here i<[{K) differs from Eq. (68) just by the exponential pre- factor induced by the scat- 
terings (cf. (103)), which yields 

Ki{K) = -2 e-i(^(^)+tH (ei^W - l) . (109) 

The expression (107) for k^{K), which is valid for 3> 1, has to approach Eq. (95) for 
K <^1. An interpolating function is given by 

Zb[K)K j\ZB) V / 

Eq. (69) is recovered for j — 1 and S — The sum — -\- is shown in Fig. 9 for 
Mh/Mi — 10"^. The agreement with the numerical result is very good. Including washout 
yields a description which correctly interpolates between the weak washout regime, rhi <^ 
m* and the strong washout regime, fhi ^ m*. 



3.3 Lower bounds on Mi 

The results for the efficiency factor are easily translated into theoretical predictions for 
the observed baryon-to- photon ratio using the relation (10). The theoretical prediction 
has to be compared with the results from WMAP [20] combined with the Sloan Digital 
Sky Survey [21], Qbh^ = 0.023 ± 0.001, corresponding to 

rj^^'' ^{6.3 ±0.3) xlO-^\ (111) 

The comparison yields the required CP asymmetry in terms of the baryon-to-photon ratio 
and the efficiency factor Hf (cf. (10)), 

^CMB / ^CMB \ 

The CP asymmetry Si can be written as product of a maximal asymmetry and an 
effective leptogenesis phase 5l [22], 

£i = £5"^ sin 5l . (113) 



29 



The connection between the leptogenesis phase and other CP violating observables is an 
important topic of current research [23]. The maximal CP asymmetry depends in 
general on Mi, rhi and, via the light neutrino masses m^, on the absolute neutrino mass 
scale m [10]. For given light neutrino masses, i.e. fixed mi and ma, Si is maximized in 
the limit mi/fhi 0, for which one obtains [24], 

£--(Mi,m) = -|-^(m3-mi) . (114) 

This expression reaches its maximum for fully hierarchical neutrinos, corresponding to 
mi = and mg = matm = A/Amf^. 

Neutrino oscillation experiments give for atmospheric neutrinos [25, 26] 

Aml^ = (2.6 ± 0.4) x lO-^eV^ , (115) 

and for solar neutrinos [27] 

Am,^„,~(7.ll^:^xlO-^)eV% (116) 

implying 

matm = (0.051 ±0.004) eV. (117) 

Thus, apart the small experimental error, matm is a fixed parameter and the resulting 
maximal asymmetry depends uniquely on Mi [24], 

^^'^-16^^;^-^^ ^To^^j lo:05^J ■ ^^^^^ 

The maximal CP asymmetry, or equivalently the maximal leptogenesis phase, cor- 
responds to a maximal baryon asymmetry 77^^. The CMB constraint 77^^ > r]^^^, 
together with Eq. (118), then yields a lower bound on the heavy neutrino mass Mi, 

Ml > Mf = 

d 3 matm «f 



^6X lO-^V V ^atm 

Note that the bound depends on the combination Tyg^^/matm whose error, after the 
WMAP result, receives a similar contribution both from r)^^^ and matm such that 

Mf^^(mi) = (6.6 ±0.8) x 10* GeV fi;f-i(mi) >4 x 10* GeV «:f-^(mi) , (120) 

with the last inequality indicating the 3a lower bound. For values of Mi <^ 10^^ GeV 
[10] one can use the results of this section for the efficiency factor, neglecting the AL = 2 
washout term. Eq. (119) then provides a lower bound on Mi which depends only on mi. 
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Figure 10: Analytical lower bounds on Mi (circles) and T; (dotted line) for mi = 0, 
^CMB = g X 10~^'^ and matm = 0.05 eV. The analytical results are compared with the 
numerical ones (solid lines). The vertical dashed lines indicate the range {msoi,matm) ■ 
The gray triangle at large Mi and large rhi is excluded by theoretical consistency (cf. ap- 
pendix A). 



Fig. 10 shows the analytical results for Mf "^(mi), based on Eq. (107) for thermal initial 
abundance (thin lines) and the sum of Eqs. (109) and (110) for zero initial abundance 
(thick lines). For comparison also the numerical results (solid lines) are shown. The 
absolute minimum for Mi is obtained for thermal initial abundance in the limit rhi 0, 
for which Kf = 1. The corresponding lower bound on Mi can be read off from Eq. (120) 
and at 3 0" one finds 

Ml > 4 X 10^ GeV . (121) 

This result is in agreement with [10] and also with the recent calculation [12]. Note that the 
lower bound on Mi becomes much more stringent in the case of only two heavy Majorana 
neutrinos [28]. The bound for thermal initial abundance is model independent. However, 
it relies on some unspecified mechanism which thermalizes the heavy neutrinos Ni before 
the temperature drops considerably below Mi. Further, the case fhi <C 10~^eV is rather 
artificial within neutrino mass models, and in this regime a pre-existing asymmetry would 
not be washed out [2]. 



31 



For zero initial abundance the lower bound is obtained for KQ{m-i) — Kpeak — 0.18, 
corresponding to rh^^ ~ 10~^eV. In this case one obtains from Eq. (120) [10], 

Mi>2xl0^GeV. (122) 

Particularly interesting is the lower bound on Mi in the favored neutrino mass range 
"^soi ^1 "^atm- This range lies in the strong washout regime where the simple power 
law scaling (108) holds. One thus obtains from Eq. (120) 

M--(mi) ~ (3.3 ± 0.4) x 10^° GeV f ^qS^) ' (^23) 



which at 3 cr implies 



Ml > 2 X 10^° GeV (jF^) ' ^ ^""^ ' ^^^^^ 

where the last range corresponds to values mgoi ^ '^i ^ ?^atm- Note that these bounds 
are fully consistent with neglecting the AL = 2 washout term, which is justified for 
Mf"" < 10^4 GeV. 

In the case of near mass degeneracy between the hghtest and the next-to-hghtest 
heavy neutrino N2, y — (M2 — Mi) /Mi <S 1, the upper bound on the CP asymmetry Si 
is enhanced by a factor ^(y) ~ l/(3y) [29, 30]. The CP asymmetry of N2 can also be 
maximal. Since the number of decaying neutrinos is essentially doubled, one obtains for 
the reduced lower bound on Mi (y < 1), 

Ml > Mr (mi, y) = , (125) 

where M™™ is given by Eq. (119). It has been suggested that for extreme degeneracies 
a resonant regime [11] is reached where e"^^^ = 0{1). In this case there is practically no 
lower bound on Mi from leptogenesis. 



3.4 Lower bound on Ti 

It is usually assumed that the lower bound on the initial temperature 7] roughly coincides 
with M™™, the lower bound on the heavy neutrino mass Mi. Here 7] can be thought of 
as the temperature after reheating, below which the universe is radiation dominated [31]. 
However, in the following we will show that this is only the case in the weak washout 
regime, i.e. for mi < m* ~ 10~^eV, whereas in the more interesting strong washout 
regime Ti can be about one order of magnitude smaller than Mi. 
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In general, the maximal baryon asymmetry is a function of both, Mi and Ti = Mi/ z\, 

r/^''"(Mi,mi,m,Zi) = ^^^"^"(Mi, mi, m) Kf(mi, ^i) . (126) 

In a rigorous procedure one would have to treat Mi and T; as independent variables and 
to determine the values zf^^ = Mi/T™™ as well as M™™ for which the CMB constraint 
^min ^ ■q'^^^ is satisfied. This will yield a value M™™|^ somewhat larger than M™™|^ 
For simplicity, we shall use the approximation M™°|^ ~ M™™|^ _q in the following. We 
then define the value z™*^^, and the corresponding temperature T™™ = M™™/^™'^^, by 
requiring that the final asymmetry rj]^^^ agrees with observation within la relative error 
of the quantity Vj^^^^ /m^i^ which controls M™™. 

In the weak washout regime, i.e. rhi < m*, one has z^'^^ ~ 1. At temperatures smaller 
than Ml, the predicted asymmetry rapidly decreases. Either, there is not enough time 
to produce neutrinos (for zero initial abundance) or the thermal abundance is Boltzmann 
suppressed (for thermal initial abundance). As Fig. 11 illustrates, for rhi < the final 
efficiency factors for zi = \ and zi \ differ by only 10%. Hence, in the weak washout 
regime one has z^^^ ~ zb — 1- 

In the strong washout regime the baryon asymmetry is predominantly produced 
around z^. The value of z^^'^ is thus given by z^^^ ~ z^ — 1.3 cr^, where is the 
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Figure 12: The function exp {—ipsiz, oo)) (cf. Eq. (106)) for mi = 10 ^eV (full line) and 
rhi = 10~^ eV (dashed line). 

width of the Gaussian which approximates the function exp {—ipsi^)) ("^f- (106)) peaked 
sX z = z-Q and whose integral between zf^^^ and infinity gives the final efficiency factor 
minus the small error that is tolerated. In Fig. 12 the function exp {—'ip^si^)) shown 
for rhi = 10~^ eV and fhi = 10~^ eV, respectively. The width of the peak is given by 
~ {iP"{zb))^^^'^ — 1.5 for fhi m*. 

One can easily write down an approximate expression for zf^^^{K) which interpolates 
between the two regimes of weak and strong washout, 



(127) 



The importance of the quantity zb becomes apparent by comparing Figs. 5 and 11. For 
instance, for fhi = matm — 0.05 eV one has zb — 8 and thus zf^^^ ~ 6, while for rhi = 
nisoi — 0.008 eV one has zb — 6 and thus z^^^ ~ 4. Clearly, for values z^ > z^^'^^K) the 
suppression of the efficiency factor becomes significant. 

From Eq. (127) one immediately obtains a lower bound on the initial temperature Ti, 



ZB-2 ■ 



:i28) 



The result is shown in Fig. 10. For small fhi < one has z^'^^ ~ 1 and consequently 
rpmm ^ ]Vf™™. Hence, in particular, the 3 a bounds (121) and (122) apply also to T;. More 
interestingly, in the favored region msoi ^ "^i ^ fna.tm (dashed box) the 3 a bound (124) 
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gets relaxed by a factor 4 to 6 and thus 

j^min = (4 X 10^ ^ 2 X 10^°) GeV (129) 

Therefore, in the favored region of fhi (dashed box), T""°(mi) is only one order of mag- 
nitude higher than the absolute minimum T; ~ 3 x 10^ GeV at rhi ~ (zero initial 
abTindance) and less than two orders of magnitude higher than the asymptotical minimum 
for nil ^ iTi-* (thermal initial abundance). This is important in view of the 'gravitino 
problem' which yields an upper bound on Ti for some supersymmetric extensions of the 
standard model. 

Comparing our results with those of [12], where the additional B — L asymmetry has 
been calculated which is produced during the reheating period at temperatures above T^^h 
and below some maximal temperature Tmax, "we find the same amount of relaxation of the 
bound on Ti = T^ch- This indicates that the relaxation is a consequence of Tb < M™™ in 
the case of strong washout rather than the existence of a non radiation dominated regime 
above Treh- 



4 Upper bound on the light neutrino masses 

We now want to study the effect of the contribution AW to the total washout. This term 
originates from the AL — 2 processes (f)l 0/ and (fxf) II with the heavy neutrinos A^i, 
N2 and A^3 in s- and t-channel, respectively. Al^ is the only term in the kinetic equations 
which is not proportional to mi but instead to the heavy neutrino mass Mi. 
At low temperatures the washout term AW reads, 

where m is the absolute neutrino mass scale, and the dimensionless constant uj is given 

by 

9x/5M^10-«GeV^ 
uj = — r-J77T ■;=—A — — 0.186 . (131) 



AW is compared in Fig. 13 of appendix A with the total washout term 

W{z) = Wo{z]mi) + AW{z]Mirff) ■ (132) 

As discussed in the appendix, there is a sharp transition to a low temperature regime 
where Al^ dominates over Wq- This transition occurs for a value ^ 1, which is 
determined by 

Wq{za) = AW{z^) . (133) 
35 



Prom Eqs. (105) and (130) one easily obtains, 

(io4v)(i)^- (-) 

In the case za ^ zb, the values of zb and of the efficiency factor at z ~ are not 
affected by AW. Since for ^ > 2;b no asymmetry is produced, the total efficiency factor 
is simply given by 

Rf{mi,Minf) = Kf{fhi) e~ -^"^^ dzAW{z) ^ (135) 

where the second factor describes the modification due to the presence of Al^. Note that 
R{ depends on rhi also via zb- For -z^ > 3 we can use the low temperature limit (130) for 
Al^, which yields 

Kf(mi, Mifn^) = Kf(mi) e"^^ (10^1^ oev) (w) . (136) 

Given the solar and atmospheric mass squared differences and a neutrino mass pattern, 
i.e. — m| > m| — mf or m| — ml < m| — ml, the dependence of m^ on mi is fixed. 
In Ref. [2] the absolute neutrino mass scale fn was used as variable. In the following we 
prefer to use instead the lightest neutrino mass mi. In the case of normal hierarchy, with 
m| — m| = Am^^jjj and m| — m^ — Amg^,, one has 

mi = ml + Aml^ + Aml,, (137) 
m^ = m^ + Am^^i, (138) 
m^ = 3m? + Am2,^ + 2Am^„i. (139) 

In the case of inverted hierarchy analogous relations hold. 
Consider now the maximal baryon asymmetry (cf. (10)), 

VB^ifhi, Ml, mi) ~ d £™^(mi. Mi, mi) Ri(mi, Mi, mi) . (140) 

In the case mi = the maximal CP asymmetry was depending only on Mi (cf. Eq. 
(118)). If mi > this is suppressed by a function f3{mi,mi) < 1 depending both on mi 
and on mi [24, 2] such that 

£r'(Mi,mi,mi) = £r'(A/i)/3(mi,mi) . (141) 

The maximal value (3 — 1 is obtained in the case mi = 0. The function (3 is conveniently 
factorized, 

/3(mi, mi) = /?max(mi) /(mi, mi) . (142) 
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The first factor, 

/3max(mi) = = ■ , (143) 

^atm [1713+ mi) 

is the maximal value of P for fixed mi, which decreases oc 1/mi for mi ^ matm- The 
function / contains the dependence on mi. 



(ms - mi + 

f(mi, mi) ^ ^- ^ . (144) 

ms — mi 

This expression describes correctly the behavior of the maximal CP asymmetry in the 
limits mi and mi oo. However, it has recently been pointed out [32] that 
Eq. (144) underestimates the maximal CP asymmetry in particular in the regime of quasi- 
degenerate neutrinos'^. For simplicity, we shall first calculate the neutrino mass bound 
using Eq. (144) and then discuss the correction. 

Let us now calculate the value Mi that maximizes rj^'^^. In the (mi, Mi)-plane this 
defines a trajectory r7^^^(mi,mi) along which r]'^^^ is maximal with respect to Mi. The 
corresponding condition. 



, (145) 



- f \ ^^V = l, (146) 



yields the relation, 

zb VlO^°GeVy VeV 

where the quantity zb is a function of mi. It is now easy to see that the ratio rj^^^/rj^^^ , 
maximized with respect to Mi, can be expressed in the following form, 

VT^'i'i'i- '"i) ^. ^, ^ ^b(wi) fjmi. nil) /Vf(mi) 

CMB — X S / I \ — 2 ' V '/ 

r]^^^^" (mi + m3j m 

where x is the constant 

X = |^— c^l.6eV3 . (148) 
Qeuj mo 

and mo = (167r/3) {v'^ / IQ^^ GeV) ~ 0.051 eV. The parameter ^ is the product 



^The expression (144) has been obtained using Eq. (22) in Ref. [2] and assuming xz = 0, which is valid 
only in the Umit mi — > 0. Note, however, that Eq. (144) approximates the maximal CP asymmetry within 
about 20% also in the quasi-degenerate case for the relevant values of mi. For quasi-degenerate neutrinos, 
with mi ~ ms < mi, one easily finds that the maximal CP asymmetry is reached for X3 ~ m\/{2m\). 
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It accounts for various factors affecting the ratio if^^ jif^^^: (1) the maximal CP asym- 
metry, (2) the atmospheric neutrino mass scale, ^atm = ?7?.atm/(0.05 eV); (3) the ob- 
served baryon asymmetry, = r{%^^ j (6 x 10~^°); (4) the variation of the strength of 
the AL = 2 washout term, and (5) the variation of the efficiency factor at small M^Wt^ . 
This parametrization of the maximal asymmetry is useful to study the dependence of the 
neutrino mass bound on the various parameters involved. 

In order to determine the absolute maximum of the asymmetry rf^^ we also have to 
find the extremum with respect to rhx and, finally, the maximum with respect to mi or, 
equivalently, the absolute neutrino mass scale. Comparison with the observed asymmetry 
jjCMB j^YiQYL yields the leptogenesis neutrino mass bound. Anticipating again that the 
maximum falls in the region of large mi, we can use the analytical expression (107) for 
Ki in the strong washout region. Since j{zB ^ 1) ~ 1, one has 

ZBirhi) Kf{rhi) ^ ^ = ■ (150) 
K mi 

Further, for large mi the function /(mi, mi) can be approximated by 

(151) 

z mj 

With this simplified expression it is easy to see that the peak is reached for 



^peak ^Jlrni (mi + mg) , (152) 



corresponding to f{'rn^^^^ mi) = 2/3. The peak value of the asymmetry is given by 

?7^*"^'^(mi) ^ 2^/^ Xfn*i 

„CMB — 03/2 Z \ 77J^ 1/2 — 2 ' 

'IB (mi + msj'^'^ mj^ m 

Anticipating (m^^^)^ ^ "^atm; one has to zeroth order in {rriatm/'iTii)^, 



(153) 



ml — —1= ~ mi . (154) 
v3 

Imposing now the CMB constraint rjB > fj^^^ we find the leptogenesis bound on the 
absolute neutrino mass scale (cf. [33, 2]), 

m,<ml,^^ae/'eV , (155) 

with 

-io»(^^^^)'"(^)"^o.m. « 
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In this last equation we used the fact that in a standard thermal history the dilution 
factor, contained in d, is given by N^^'^/N* — g*/grec with g^ec = 43/11 the number of the 
(entropy) degrees of freedom at recombination. Combining Eqs. (154) and (152) one finds 



m 



peaJj 



~m^,,k-0.17eV, (157) 



which is consistent with the approximation of strong washout used in Eq. (150). From 
Fig. 5 one then reads off z^^"^^ ~ 10. Together with (146) this yields for the peak value of 
Ml, 

/ ~ peak \ 

Mf^ = "^^^^1 ^{ ~ 2 X 10^^ r'/' GeV . (158) 

It is straightforward to go beyond the zeroth order in {matm/'miy . In the case of 
normal hierarchy the lightest neutrino mass bound is given by 

<.t^<.{^-'i,J^\, (159) 
whereas in the case of inverted hierarchy one has 

-l,inv - -,eak ' ) , (160) 

which yields m^^^^ — "^i^nt — "^peak ~ 0.005 eV. In order to obtain numerical results for 
the upper bounds on the light neutrino masses one has to specify the baryon asymmetry 
and the neutrino mass squared differences. For rj^^^ we use the WMAP plus SDSS result 
(111), while the value for matm is given by the Eq. (117). Since mpeak oc {ml^^/rjcMBY^'^, 
the experimental error on m^^^^ is about 5%. Setting all other parameters — 1, one 
finds for the central value ^ = Catm/C?? — 0-95. We then obtain 

The corresponding 3a upper bounds on the neutrino masses are for normal hierarchy, 

mi, m2 < 0.13 eV , ma < 0.14 eV , (162) 

and correspondingly in the case of inverted hierarchy, 

mi<0.13eV, ms, mg < 0.14eV . (163) 

These analytical bounds are consistent with the numerical results obtained in [2], if one 
accounts for the different parameters, matm = 0.05 eV and tjcmb — 3.6 X 10-1° (C- 1-7), 
and the over-estimate of the washout term Wo by 50%. 

39 



The bound on the maximal CP asymmetry derived in [32] corresponds in the relevant 
range of large mi and rhi to the function (cf. (144)), 



Repeating the above analysis one finds that the peak of the asymmetry is shifted to 
-peak ^ y2mi, with /P'^^'^ ~ l/V^. From Eqs. (150), (154) and (157) one then con- 
cludes that the neutrino mass bound is relaxed by the factor (3^/^/4)-^/^ ~ 1.07, i.e. 7%, 
corresponding to an increase of the neutrino mass bound by 0.01 eV . 

An important correction arises from the dependence of the neutrino masses on the 
renormalization scale /i. The only low energy quantity upon which m\^^^ depends is the 
atmospheric neutrino mass scale matm- Hence, there arc two competing effects: the first 
one is the running of matm from the Fermi scale /i = uiz to the high scale /i ~ Mi 
(~ 10^^ GeV), the second one is the running of m\^^^ from ~ Mi down to /i = raz- In 
the standard model the light neutrino masses scale uniformly under the renormalization 
group. Since mpeak oc y^m^^m, the first effect then gives a correction that is only half of the 
second one. Renormalization group effects make the bound more restrictive [34] . In order 
to have an upper bound, we have to choose those values of the parameters that produce 
the smallest effect. This corresponds to choosing the lowest Higgs mass compatible with 
positive Higgs self-coupling at ~ 10^^ GeV, which is about 150 GeV. The atmospheric 
neutrino mass scale is then increased by about 40% [34] and the bound gets 20% weaker 
at // ~ Ml, but 20% more restrictive at // = m^. Gombining the effect of radiative 
corrections and the larger CP asymmetry (Eq. (164)), we finally obtain from Eq. (162) 
and (163) at 3(7, 



which, thanks to cancellations among different corrections, agrees with the bound obtained 



It is important to realize, however, that there are corrections of the same order as 
those discussed above which cannot be treated within the present framework. It is usu- 
ally assumed that in leptogenesis first a lepton asymmetry is produced, which is then 
partially transformed into a baryon asymmetry by sphaleron processes. However, this 
picture is incorrect [35] . The duration of leptogenesis is about two orders of magnitude 
larger than the inverse Hubble parameter when it starts. Since many processes in the 
plasma, in particular the sphaleron processes, are much faster, the generated asymmetry 
is 'instantaneously' distributed among quarks and leptons. Hence, the chemical potentials 
of quarks and leptons are changed already during the process of leptogenesis. A complete 




(164) 



rrii < 0.12 eV , 



(165) 



in [2]. 
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analysis has to take into account how the contributing 'spectator processes', which are in 
thermal equilibrium, change with decreasing temperature (cf., e.g., [36]). In [35] is has 
been estimated that spectator processes reduce the generated baryon asymmetry by about 
a factor of two. Hence, there is presently a theoretical uncertainty of at least 0.02 eV on 
the neutrino mass bound. 

Finally, it has to be kept in mind that our whole analysis is based on the simplest 
version of the seesaw mechanism with hierarchical heavy Majorana neutrinos. The lepto- 
genesis neutrino mass bound can be relaxed if the heavy Majorana neutrinos are, at least 
partially, quasi-degenerate in mass. In this case the CP asymmetry can be much larger 
[29, 30] than the upper bound used in the above discussion. This possibility has to be 
discussed in the context of realistic models of neutrino masses. Further, if Higgs triplets 
contribute significantly to neutrino masses the connection between baryon asymmetry 
and neutrino masses disappears entirely. Different relations between neutrino masses and 
baryon asymmetry are also obtained in non-thermal leptogenesis [37, 22]. 

5 Towards the theory of thermal leptogenesis 

The goal of leptogenesis is the prediction of the baryon asymmetry, given neutrino masses 
and mixings. The consistency of present calculations with observations is impressive, but 
so far it is not possible to quote a rigorous theoretical error on the predicted asymmetry, 
which is a necessary requirement for a 'theory of leptogenesis'. 

The generation of a baryon asymmetry is a non-equilibrium process which is generally 
treated by means of Boltzmann equations. This procedure has a basic conceptual problem: 
the Boltzmann equations are classical equations for the time evolution of phase space 
distribution functions; the involved collision terms, however, are zero temperature S- 
matrix elements which involve quantum interferences in a crucial manner. Clearly, a full 
quantum mechanical treatment is necessary to understand the range of validity of the 
Boltzmann equations and to determine the size of corrections. 

A first step in this direction has been made in [38] where a perturbative solution of the 
exact Kadanoff-Baym equations has been constructed. To zeroth order, for non-relativistic 
heavy neutrinos, the non-equilibrium Green functions have been obtained in terms of 
distribution functions satisfying Boltzmann equations. Note that in the favoured strong 
washout regime the decaying heavy neutrinos are indeed non-relativistic. It is instructive 
to recall the various corrections. There are off-shell contributions, 'memory effects' related 
to the derivative expansion of the Wigner transforms, relativistic corrections and higher- 
order loop corrections. All these correction terms are known explicitly, but their size 
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during the process of baryogenesis and, in particular, their effect on the final baryon 
asymmetry have not yet been worked out. 

Recently, thermal corrections have been studied [12]. They correspond to loop cor- 
rections involving gauge bosons and the top quark. At large temperatures, T > Mi, the 
processes in the plasma and the CP asymmetries change significantly if thermal masses 
are treated as kinematical masses in the evaluation of scattering matrix elements [12]. On 
the contrary, thermal corrections are small if they are only included as propagator effects 
[14] . To clarify this issue is of general importance for the treatment of non-equilibrium 
processes at high temperatures. 

The effect of all these corrections on the final baryon asymmetry depends crucially on 
the value of the neutrino masses. Large thermal corrections would modify the asymmetry 
at temperatures above Mi. This affects the final baryon asymmetry only in the case of 
small washout, i.e. rhi < m*. In the strong washout regime, rhi > m*, which appears to 
be favored by the current evidence for neutrino masses, the baryon asymmetry is generated 
at a temperature Tg < Mi. In this case thermal corrections are small. Correspondingly, 
the recently obtained bounds on light and heavy neutrino masses [2, 12, 32] are all very 
similar. 

The final value of the baryon asymmetry is significantly affected by 'spectator pro- 
cesses' [35] which cannot be treated based on the simple Boltzmann equations discussed 
in this paper. It has been estimated that this effect changes the baryon asymmetry by a 
factor of about two, leading to a theoretical uncertainty of the leptogenesis neutrino mass 
bound of about 0.02 eV. Clearly, to obtain a more accurate prediction for the baryon 
asymmetry requires a considerable increase in the complexity of the calculation. 

An important step towards the theory of leptogenesis would be a systematic evalua- 
tion of all corrections to the simple Boltzmann equations in the 'easy regime' of strong 
washout where fhi > m*. One could then see where this approach breaks down as rhi 
decreases and Tb approaches Mi. On the experimental side, information on the absolute 
neutrino mass scale m, and therefore on rhi > "^i, is of crucial importance. Maybe, we 
are lucky, and nature has chosen neutrino masses in the strong washout regime where 
leptogenesis works best. 
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Appendix A 



A crucial and delicate point in setting up the Boltzmann equations for leptogenesis is the 
subtraction of the real intermediate state contribution (RIS) from the 2 — > 2 scattering 
processes [4]. Without this subtraction, decays and inverse decays lead to the generation 
of a lepton asymmetry in thermal equilibrium, in contradiction with general theorems. 

In order to explicitly split the 2 — > 2 scattering processes into RIS and remainder 
one has to calculate the AL = 2 processes, including one-loop self-energy and vertex 
corrections in the resonance region. This calculation has been carried out in [39] where 
the relevant results are given in Eqs. (68) - (85). 

Let us consider for simplicity the case s, Mf -C Ml,Ml (s, t and u are the usual 
Mandelstam variables). For our purposes it is sufficient to study the averaged matrix 
element squared, 

\Mm^mi,= f du\Mm^m\ {im) 

J —s 

where the integral over u corresponds to the integral over the final state lepton angle, i.e. 
a partial phase space integration. 

In order to study the resonance region the diagonal part of the self-energy was re- 
summed in [39] whereas the off-diagonal part was treated as perturbation in the Yukawa 
couplings h. In the 2^2 scattering amplitude the free propagator is then replaced by a 
Breit-Wigner propagator. 



s - Mf + iMiFi Ml \Di{x) Ri{x) J ' 

where 

1 x-1 1 



Di{x) (a;-l)2 + c2' Ri{x) {x - If + 



(168) 



and 



with 



D^{xY Ri{xf cRi{x) ' 

s _ Ti _ 1 

Ml ' ^~ Ml ~ 87 

here hkj is the Yukawa coupling of Nj to lk4> 



1111 , , 

+ -fTT^ = -TwT ' (169) 
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The averaged matrix elements are then given by the following expression [39], 

|A^([^^Z0)|L = 2s'Y,lA^-B,^-Qj-isY,{Dijk + E,^k)], (171) 

ij \ k / 

\M(l(l>^ml = 2s^J^lAj + Bij + Cij + 4sJ2{Dijk + Eijk)], (172) 

ij \ k / 

where we have only shown terms contributing to the subtraction of the RIS part as well 
as the leading order off-shell part, e.g. CP conserving one-loop corrections have not been 
included. Aij and Bij represent the various Ni-Nj s-channel interference terms. Up to 
terms 0{h^) they are — 2, 3) 



.2 W 1 1 



+ = -2 Re{i.?J^i- , (174) 

Bu = B,j + B,, = 0, (176) 

B.. + B., = 2 Im{i.,^J^i- ; (177) 

the Ni-Ni s-channel terms with self-energy and vertex corrections, respectively, read 
(A; = 2, 3) 



1 1 



= -2x lv^{Kl,}^^ ( 7^ + 151 ) , (178) 



where 



m = Vy(^i-(i + y)in(^i±^) 



^ +o(-] ; (180) 



2^/y \y 
finally, the s-u-channel interference term is 

Eiik — —Duk ■ (181) 
From these equations one reads off 

i^m + ^m = 0, (182) 
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and, for s — , 

Bu + Bn + Cu^O . (183) 

Hence, the CP asymmetry of the full 2 — > 2 cross section vanishes to 0{h'^). The 'pole 
terms', corresponding to Ni-Ni s-channel contributions, are cancelled by on-/off-shell 
s-channel interferences (self-energy) and s-channel/u-channel interference (vertex correc- 
tion). Off-shell, the corresponding cancellations take place to 0{h^) [40]. 

As an unstable particle the heavy neutrino A^i is defined as pole in the 2 — > 2 scattering 
amplitude 



{Ni\U) . (184) 



The residue yields the decay amplitude and, in particular, the CP asymmetry. The RIS 
term can then be identified as the squared matrix element in the zero- width limit, 

\M{l4>^mlis = hm 2^2 ||An-X^(Cifc + 4sDnfe)j 

= IGTr^XiiMi' (1 + 2£f + 2eX) 5{x - 1) , (185) 

where and e\ arc the familiar CP asymmetries due to mixing and vertex correction, 
respectively (Mi < M^,M^), 

' 8x ^ ifii ' ^ ' 

k 

y _ 1 ^ l^{Kl,} (Ml\ 

In [4] it has been shown that the subtraction of the RIS term, which corresponds to 
the replacement {l/Df + 1/Rf) {l/Dj + l/Rj) - {n/c)6{x - 1), leads to Bohzmann 
equations with the expected properties, which have the equilibrium solution A^i — N^, 
Nb-l = 0. 

It is now straightforward to write down the subtracted matrix element squared. Keep- 
ing only terms 0{h'^), where the zero- width limit can be taken for B^ and Enk, one 
obtains the simple expressions = 2,3), 

iMiU^ml, = \M^l=2\1 + \Mal=2\'., (188) 
\M{l(l>^m\l, = \Mal=2\1 - \Mal=2\1 , (189) 
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with 



\M 



AL=2\ 



+ 



2s' 



-"■11 



2 2 

+ 2 



1 + 1 ln(x + 1) + 



-6 yRe(x2 ) 1 



MiMi 



1 1 
- + 



X 2Di 



xDi 

(^+ l)i / .^ 

- — ^In x + 1 

X'^ 



+ 3VRe(ii'J)-i-| , 
^ ^ MiMi I 



(190) 



AL=2 _ 



= -327r'XnM2(£f + £]")5(x - 1) . 



(191) 



Note that the subtracted squared matrix element, contrary to the unsubtracted one, 
violates CP. To leading order in the coupling this part contributes only on-shell, and 
it is Oih?) suppressed with respect to the leading Born term. Away from the pole, for 
s <S Ml, one has 



6.^^Re{i^5}-^ 



MM, 



tr(mj,m,,) , 



(192) 



which is the crucial term leading to the upper bound on the neutrino masses [10]. 

We also have to take into account the AL = 2 process II — > 4>4>- The corresponding 
matrix elements reads 



\M 



AL=2,t 



2s' 



Ml 



X + 1 x{x + 2) 



ln(x + 1) 



+6 Re iKl) — ln(a; + 1) + 3 V Re (kD — ^ I . (193) 



For small center of mass energies one again obtains 



\M.AL=2,t\ = -^tr{mlm„) . 



(194) 
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For the derivation of the upper bound on the hght neutrino masses one needs the 
maximal CP asymmetry for given Mi, mi and fn. In this case also the complete AL = 2 
matrix element depends just on these three variables. This is easily seen in the flavor 
basis where the Yukawa matrix h connects light and heavy neutrino mass eigenstates. 
The matrix 



\/rniM^ 



■A 



ij 1 



is then orthogonal, = Vt^Q, — I [41], which implies 



(195) 



(196) 



Using Re{i^ij} = Re{ir?i}, this imphes for the interference term appearing in Eq. (190), 



K 



11 



Ml 



The conditions 



Re{hl^} = Re{/i^i} = , Re{/i?i} 



miMi 



(197) 



(198) 



yield a good approximation for the maximal CP asymmetry [2]. The difference to the 
maximal CP asymmetry [32] can then be treated as a perturbation, as discussed in sect. 4. 
Eq. (198) then implies for the interference term 



(199) 



Inserting this expression in Eq. (190) one obtains for the AL = 2 matrix element in the 
case of maximal CP asymmetry. 



|MaL=2U = ^ 



mi 



-"2 + -^2 - -8{X - 1) + - 1 + — 



X 



(200) 



x + l 



+6 (mi — mi) 



x + l 1 

h 



ln(a; + 1) 
(x+l) 



X 



2Di 



x^ 



ln{x + 1) 



+ 3m2 



For the process II — > (fxf) one obtains in the case of maximal CP asymmetry. 



\M 



AL=2,t\ 



2£2 



mi 



x + l 



3 + 



x{x + 2) 



+6 (mi — mi) 
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1 ln(a; + 1) 

X 



ln{x + 1) 



+ 3m^ 



(201) 



For small energies, s <^ Ml, these matrix elements again reduce to 

6s2 



\Mal=2\1 = |MAL=2,t|' = , (202) 

whereas for intermediate energies -C s <^ -^23 finds 

|A<AL=2|+ = |A<AL=2,t|' = — (m? + - 2m2) . (203) 



Following [4, 5], it has been standard practice [6]-[ll] to determine |A^|s„b by comput- 
ing the Born diagrams for the 2 — > 2 process with Breit-Wigner propagator and dropping 
the imaginary part of the propagator, since in the zero width limit 

1 1 



S-MI + iMiFi Ml \ Di {x) Ri {x) , 

"^-^ -i7rS{s -Ml). (204) 

Recently, it has been pointed out that this procedure is not correct [12]. In a toy model, 
the same conclusion has been reached in [42]. Indeed, the described procedure leads to a 
subtracted squared matrix element which contains terms oc 1/Df [5], implying — 
0{h?) on-shell, in contradiction with Eq. (190). The zero-width hmits of the squared 
amplitude and the squared imaginary part are different. This was overlooked in the past, 
leading to an overestimate of the washout rate due to inverse decays by 50% [12]. The 
RIS term has to be subtracted from the full 2 — > 2 cross section, not just from the Born 
cross section, in order to obtain the crucial CP violating contribution proportional to £1^. 

Let us now consider the Boltzmann equation for the density of lepton doublets, as- 
suming kinetic equilibrium, 

^ + mni = "^^-^^N, ^ lct>) - ^^^\lct> ^ N,) (205) 
at n^^ rii 

"■I "■I 

+ 7*^^(00 ^ U) - 7'^'^(^^ ^ #) ■ 

^Note that in Ref. [12] -f'}^~^^''^\LH LH) is CP violating, using a given CP asymmetry ejvi not 
determined by the 2 — > 2 processes. It is different from the on-shell part of {LH — > LH) which, as a 

tree-level rate, conserves CP. As discussed above, the correct subtraction term is obtained from the full 
reaction rate including vertex and self-energy corrections to 0{h'^) by separating the on-shell part from 
the interference terms. This procedure automatically yields the correct CP asymmetry e^i- 
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Here, j^"^ are the usual reaction densities in thermal equilibrium and we have assumed 
that the Higgs doublets are in thermal equilibrium, neglecting their chemical potential. 
The CP asymmetry £i is defined in such a way that 

t'^W ^ ^0) = ^ A^i) = ^^7iv, , (206) 

^ ^"^) = ^ ^i) = ^^7iv, . (207) 

Further, for the 2 — > 2 processes we have 

iZiU^m = 71W - ^^m, , (208) 
7:i(/0^/0) = 7lW + ^^i^^i' (209) 

Introducing a lepton, or B — L asymmetry, 

ni = nl"^ - ^^B-L , nj = n)"' + ^^b-l , (211) 

assuming ub-l — nj— rii — 0{ei), and keeping only terms 0{ei), one obtains the kinetic 
equation for the B — L asymmetry 



dnB-L , „ „ f ^\ riB-L / 1 , , (^^^. 
h SHub-l = -Si - I J 7iVi - \2^^^ ^^^=^ ' ' ( ) 



where 

7AL=2 = 2 7l\^2,+ + 2 7rL=2,t ■ (213) 

The CP violating part of 7g^,^ yields the term +£i7Ari which guarantees that for — "rf^i 
no asymmetry is generated. Note that the old procedure for subtracting the RIS part of 
the 2^2 process would have led to a contribution ^'Jni in the washout term rather than 
^^Ni [12]. Neglecting the off-shell contribution 7al=2, using the relation 

7;vi ^rf^zH D = 2nl'' zH , (214) 

and changing variables from tio z — Mi/T and from number densities to particle numbers 
in a comoving volume (cf. [10]), one obtains the Boltzmann equation (13). 

In order to obtain the kinetic equation (212) the correct identification of the RIS term 
is essential. It is therefore of crucial importance to derive this equation from first prin- 
ciples. In the case of non-relativistic heavy neutrinos, i.e. T < Mi, this has been done 
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in [38]. Note that in the strong washout regime, where Tb < Mi, the decaying neutrinos 
are indeed non-relativistic. Eq. (49) of [38] gives the analogue of (212) for Boltzmann 
distribution functions rather than the integrated number densities. The starting point 
of this derivation are the Kadanoff-Baym equations which describe the full quantum me- 
chanical problem. Leptogenesis is then studied as a process close to thermal equilibrium. 
As a consequence, the deviations of distribution functions from equilibrium distribution 
functions, Sfiy{t,p) and 6fi{t,k) appear from the beginning^. For simplicity, in Eq. (49) 
of ref. [38] the contribution to the washout term from interferences with the heavy neu- 
trinos N2 and A^3 has been neglected. Otherwise the result is identical to Eq. (212). In 
particular, the relative size of the driving term for the asymmetry, which is proportional 
to El, and the washout term due to inverse decays agrees with (212). It is important to 
derive the Boltzmann equations and the reaction densities within a full quantum mechan- 
ical treatment also for relativistic heavy neutrinos, in particular in the resonance region 
T ~ Ml. 

It is instructive to discuss the different contributions to 7al=2, the reaction density 
corresponding to the averaged matrix element squared 



\M 



AL=2\ 



2 \MaL=2\+ + 2 \MAL=2,tf 



4s2 



1 1 



327r 



2miv 



1 



6 + - + 

X 



2 1 

x + 1 ^ 'd[ 



Ml x2 
2 



5{x - 1) + 24m2 (215) 



3 + 



X 



x^ 



6x + 3 + 



— 6 mi 



x + 2 
1 



+ 



x + 1 

(2x + 



ln(a; 1) 



\n{x 1) 



where wc have again assumed the relation (199). The different contributions to the 

washout term W are shown in Fig. 13. The term proportional to m^, as well as the 

contributions from the last three lines at high temperatures, Mi ^ T <^ M2,3, give a 

simple power law behavior, corresponding to Eqs. (202) and (203). At low temperatures, 

the term proportional to m\ rapidly approaches zero and becomes negligible. 

It can be seen clearly that for z < 30 the contributions from the first two terms in 

Eq. (215) cancel each other to a very good approximation, corresponding to the subtrac- 

^The part of the Lagrangian involving left- and right-handed neutrinos has a U(l) symmetry which 
implies dfi{t,k) = —Sf^{t,k). Here we have assumed that due to the other interactions in the standard 
model Sf^{t, k) = 0. 
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Figure 13: Absolute values of the different contributions to the washout term W from 
the averaged squared matrix element in eq. (215) for Mi = 10^° GeV, fn = 0.05 eV, 
mi = m/v^ and fhi = 0.03 eV. The solid line is the term proportional to l/Df + l/Rf, 
the black squares the contribution from the delta function, the dotted line the term 
proportional to in the second and third lines of eq. (215), the dashed line the term 
proportional to ml, and the dashed-dotted line is the term proportional to m^. 

tion of RIS contributions. However, the term proportional to 1/ + 1/R1 has a different 
low temperature limit than the delta function and cancels against the term in the second 
and third lines of Eq. (215) for z > 30. 

Finally, this discussion is only applicable if off-shell and RIS contributions can be 
separated. This is related to the usual approximation that the right handed neutrinos 
can be considered as asymptotic free states, i.e. that one can write down a Boltzmann 
equation for them, which is the case if their width is small, i.e.. 



where S is some constant smaller than one. This translates into the following condition 
for mi. 




(217) 
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We have checked numerically that the separation of on-shell and off-shell contributions 
works well, as long as Eq. (217) with 5 = 0.1 is fulfilled. 



Appendix B 

The scattering rates are expressed through the reaction densities rates and these, in turn, 
through the reduced cross sections. 



where we introduced the following integrals 



<j>,tis) 



The reduced cross sections can be written in the following form [8] : 

3 Q!^ Ml rhi 



cr<i>,t{s) 



47r 

with X = ip / z"^ , = ml/v'^ and where we defined the following functions: 

x-l 



X 



X - 2 + 2ah 1 - 2ah , f x - 1 + an 
H — In 



X -1 + Qh 



x-l 



ah 



X — 1 



X 



with ah = (Mh/Mi)^. The functions f(j,,t{s){z) are then defined as: 



z^Ko(z) 



(218) 



(219) 



(220) 



:22ii 



(222) 



(223) 



and in this way the Eq. (73) for St follows. 

Similarly to the Eq. (25) for K2{z), the modified Bessel function K3{z) can be ap- 
proximated by the analytical expression 



945 35 



For Mh/Mi = IQ-^ and small K, Zeq is well described by 



(224) 



^eq = 0.4 + 1.3 ln(l + X-o-8«) . 



(225) 
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